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SUMMARY  AMD  RECOMMENDATIONS 

The  purpose  of  chis  study  was  to  analyze  the  time  series  of  the  second 
destination  tonnage  and  to  suggest  models  to  forecast  them.  The  time  series 
models  to  be  attempted  were:  exponential  smoothing,  polynomial  regression,  and 
Box-Jenkins.  A  structural  "regression"  model  using  overseas  flying  hours  as 
the  independent  variable  was  also  to  be  attempted. 

Table  1  gives  the  assortment  of  models  we  used.  It  can  be  ascertained 
that  we  examined  more  models  that  originally  contemplated.  The  reason  was 
some  special  problems  with  the  sealift  series:  seasonality  in  some  cases;  and 
the  confounding  of  the  series  by  the  change  of  the  fiscal  year  in  1977. 

An  inspection  of  Table  1  reveals  that  a  range  of  models  can  be  used  to 
forecast  tonnage,  the  models  themselves  spanning  the  very  simple  to  the 
extremely  sophisticated.  The  simplest  model  used  was  that  of  the  sample  mean. 
This  was  then  extended  to  provide  high-low  estimates.  The  most  complicated 
model  used  was  the  Wiener- Kolmogorov  method,  an  autoregressive  spectral  analysis 
technique . 

In  recommending  one  forecasting  model  for  any  one  tonnage  series,  we  have 
to  adhere  to  the  principle  of  cost-effectiveness.  That  is,  if  we  recommend  a 
forecasting  technique  to  you,  we  wish  it  to  be  relatively  accurate  and  at  the 
same  time  we  would  like  it  to  cost  you  as  little  as  possible.  This  principle 
has  led  us  to  recommend  the  two  types  of  models  listed  in  Table  2:  single 
exponential  smoothing  for  nine  of  the  eleven  tonnage  series;  and  time  series 
decomposition  for  the  remaining  two.  There  is  no  implementation  cost  in  terms 
of  computer  programming  or  training  In  following  these  recommendations.  For 
exponential  smoothing,  one  need  only  store  two  values:  the  current  tonnage 
and  the  forecast  for  the  current  period.  As  for  time  series  decomposition. 


one  need  only  use  the  seasonal  Indices  that  we  have  already  calculated.  These 
indices  may  be  adjusted  from  time  to  time,  according  to  the  judgment  of  the 
manager. 
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TABLE  2 

FORECASTING  MODELS  RECOMMENDED 


Series 

Name 

Model 

Forecasting  Model 

# 1 

PAC/AIR 

Single  Exponential  Smoothing 

St+1  *  0,60Xt  +  °-40St 

n 

ALA/AIR 

It 

St+1  *  °-24Xt  +  °-76st 

03 

NE/AIR 

It 

St+i  -  0.65Xt  +  0.35St 

0  4 

EUR/AIR 

ft 

S  i  -  0.49Xt  +  0.51St 

05 

SOUTH/AIR 

It 

S  x  -  0.10Xt  +  0.90St 

#6 

PAC/SEA 

ft 

S  +i  -  0.99Xt  +  0.01St 

07 

ALA/ SEA 

It 

St+1  “  °'42Xt  +  °'58St 

its 

NE/SEA 

Time  Series  Decomposition 

See  Annex  F  for  the 
seasonal  index 

It9 

EUR/ SEA 

Single  Exponential  Smoothing 

S  j  =  0.42X  +  0. 58S 

t+1  t  c 

It 10 

SOUTH/SEA 

Time  Series  Decomposition 

See  Annex  F  for  the 
seasonal  index 

till 

CONUS 

Single  Exponential  Smoothing 

S  .  =  0.60X  -  0.40S 
t+i  c  c 

Symbols: 

St+1  Is  the  forecast  made  in  period  t  for  period  t  +  1 
Xfc  is  the  actual  value  realized  in  period  t 
Sfc  is  the  forecast  made  in  period  t  -  1  for  period  t 
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ANNEX  A 

PLOT  OF  THE  MODIFIED  DATA  AND  THE  DEVELOPMENT 
OF  LOW-HIGH  FORECASTS  FOR  THE  STATIONARY  SERIES 


Exhibit  A  provides  a  plot  of  the  tonnage  data  versus  time.  This  is  given 
for  the  eleven  sets  of  tonnage  series  provided. 

The  plot  of  the  data  leads  us  to  recommend  forecasting  models  to  use  in 
each  case.  Series  #1  (PAC/AIR)  appears  as  a  random  walk;  i.e.,  it  is  a 
stationary  process.  Hence  ,  a  ninety-five  percent  probability  interval  giving 
low-high  forecasts  would  be  appropriate.  We  shall  also  use  other  models  for 
forecasting  this  series:  single  exponential  smoothing;  and  some  selected  Box- 
Jenkins  models.  Similarly,  Series  II 8  (NE/SEA)  exhibits  a  seasonal  component; 
witness  the  spikes  during  the  July/August/September  quarter.  Consequently,  we 
would  recommend  the  use  of  a  decomposition  model.  And  so  on! 

The  recommendations  concerning  the  models  to  use  are  listed  in  Table  A.l. 
Also  listed  are  relevant  information  including  the  number  of  observations,  N. 

We  have  censored  the  last  three  observations  of  each  monthly  series  for  purposes 
of  comparing  forecasts  to  actuals.  These  observations  would  be  for  July, 

August,  and  September  of  1978.  For  the  quarterly  data,  we  censored  the  third 
and  fourth  quarters  of  1978. 
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Series 

PAC/AIR 


ALA/AIR 


NE/AIR 


EUR/AIR 


SOUTH/ AIR 


PAC/SEA 


ALA/SEA 


NE/SEA 


TABLE  A.l 

MONTHLY  DATA:  TIME  SERIES  MODELS 


Series 

No. 

Start  Date 

End 

Date 

N 

Visual 

Identification 

Models  to  Use 

1 

July  1976 

June 

1978 

24 

Random  Walk 

a.  Probability  Interval 

b.  Single  and  Double 

Exp.  Smoothing 

c.  Selected  Box-Jenkins 
Models 

2 

July  1972 

June 

1978 

72 

Random  Walk 

or 

Quadratic 

a.  Probability  Interval 

b.  Single,  Double,  and 
Triple  Exp.  Smoothing 

c.  Optimal  Box-Jenkins 
Model 

3 

July  1976 

June 

1978 

24 

Random  Walk 

a.  Probability  Interval 

b.  Single  Exponential 
Smoothing 

c.  Selected  Box-Jenkins 
Models 

4 

July  1972  June 

(Modified  Data) 

1978 

72 

Random  Walk 

a.  Probability  Interval 

b.  Single  and  Double 
Smoothing 

c.  Optimal  Box-Jenkins 
Model 

5 

Jan.  1974 

June 

1978 

54 

Random  Walk 

a.  Probability  Interval 

b.  Single  Exponential 
Smoothing 

c.  Optimal  Box-Jenkins 
Model 

6 

1st  Qtr . , 
1973 

(Corrected 

2nd  Qtr., 
1978 

Data) 

23 

Curvilinear 

a.  Single  &  Double 

Exp.  Smoothing 

b.  Selected  Box-Jenkins 
Models 

7 

II 

(Corrected 

II 

Data) 

23 

Curvilinear 

a.  Single,  Double,  and 
Triple  Exp.  Smoothing 

b.  Selected  Box-Jenkins 
Models 

8 

II 

(Corrected 

II 

Data) 

23 

Seasonal 

a.  Time  Series  De¬ 
composition  Without 
Trend 

TABLE  A.l  (Continued) 

MONTHLY  DATA:  TIME  SERIES  MODELS 
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Series 

Series 

No. 

Start  Date 

End  Date 

N 

Visual 

Identification 

Models  to  Use 

EUR/SEA 

9 

1st  Otr . ,  2nd  Qtr . , 

1973  1978 

(Corrected  Data) 

23 

Straight  Line 

a.  1st  Order  Polynomial 

b.  Single,  Double,  and 
Triple  Exp.  Smoothing 

c.  Selected  Box-Jenkins 
Models 

SOUTH/SEA 

10 

tf 

(Corrected 

II 

Data) 

23 

Complex  with 

Some  Seasonality 

a.  Decomposition  Model 

b.  Sinusoidal  Model 

CONUS 

11 

July  1974 

March  1978 

45 

Random  Walk 

or 

Quadratic 

a.  Probability  Interval 

b.  Single,  Double  and 
Triple  Exp.  Smoothing 

c.  Optimal  Box-Jenkins 
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ANNEX  B 

HIGH-LOW  FORECASTS  FOR  THE  RANDOM  WALK  MODEL 

We  can  make  the  assumption  that  some  of  the  tonnage  series  conform  to  the 
random  walk  model.  These  are:  Series  111  -  II 5  and  Series  #11.  The  rest  of 
the  Series,  II 6  -  It  10,  do  not  exhibit  stationarity . 

If  we  assume  that  the  data  for  the  stationary  series  are  normally  distributed, 
we  can  construct  ninety-five  percent  probability  intervals  for  the  population 
values  of  the  series.  The  lower  limit  of  the  interval  would  be  considered  as  a 
low  forecast;  the  higher  limit  of  the  interval  would  be  a  high  forecast.  Table 
B.l  gives  these  limits,  calculated  using  the  expression  X  +  1.96s.  (Note  that 
this  is  different  from  constructing  a  ninety-five  percent  confidence  interval  for 
the  population  mean,  in  which  case  we  would  use  the  expression  X  +  1.96  s//n.) 

If  we  wished  one  value,  rather  than  a  high-low  range  of  values,  the  best 
forecast  would  be  the  mean  of  the  series,  X.  We  compare  the  mean  of  the  series 
with  the  actuals  for  July  1978  in  Table  B.2.  We  see  that  the  percentage  error 
is  relatively  small.  What  is  more  important  is  that  all  the  actuals  for 
July  1978  are  within  the  low-high  limits  given  in  Table  B.l. 
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TABLE  B.l 


HIGH- 

-LOW  FORECASTS  FOR 

SELECTED  TONNAGE  SERIES 

; 

RANDOM 

WALK 

Series  No. 

Series 

Arithmetic 
Mean.  X  (tons) 

Standard 

Deviation,  s  (tons) 

Forecast 

Low 

(tons) 

High 

•  1 

PAC/AIR 

1,888.5 

178 

1,540.4 

2,236.6 

2 

ALA/AIR 

617.1 

98.5 

424.0 

810.2 

3 

NE/AIR 

273.8 

61.4 

143.4 

394.1 

4 

EUR/AIR 

2,732.8 

223.0 

2,295.1 

3,169.3 

5 

SOUTH/ AIR 

106.4 

27.0 

53.5 

159.4 

6 

PAC/SEA 

85,113.0 

31,236.0 

Not  Applicable 

7 

ALA/SEA 

4,727.2 

1,761.0 

Not  Applicable 

8 

NE/SEA 

2,946.9 

3,833.0 

Not  Applicable 

9 

EUR/SEA 

66,118.0 

13,060.0 

Not  Applicable 

10 

SOUTH/SEA 

4,674.5 

1,806.0 

Not  Applicable 

11 

CONUS 

1,897.8 

304.0 

1,302.0 

2,493.6 

18 


TABLE  B.2 


PERCENTAGE 

ERROR  IN  ONE-STEP 

AHEAD  EORECASTS 

USING 

THE  RANDOM  WALK 

MODEL 

Series  No. 

Series 

Arithmetic 
Mean.  X  ( t ons) 

Actual  for 
July  '78 

Percentage 
_ Error 

1 

PAC/AIR 

1888.5 

2,011 

6.9 

2 

ACA/AIR 

617.1 

565 

-  9.2 

3 

NE/AIR 

273.8 

228 

-20.1 

4 

EUR/AIR 

2732.8 

2,506 

-  9.1 

5 

SOUTH/AIR 

106.4 

97 

-  9.7 

11 

CONUS 

1897.8 

1,475 

CO 

Csl 
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ANNEX  C 

EXPONENTIAL  SMOOTHING  MODELS 

This  annex  presents  the  exponential  smoothing  forecasting  models  that 
performed  best  In  making  one-step-ahead  forecasts.  Table  C.l  lists  the 
eleven  tonnage  time  series  in  question,  the  Brown  exponential  smoothing  model 
used,  the  corresponding  smoothing  constant,  the  one-step-ahead  forecast,  and 
the  percentage  forecast  error.  Because  some  series  such  as  NE/SEA  (Series  #8) 
did  not  fare  well  using  Brown's  model,  another  exponential  smoothing  procedure 
which  incorporates  trend  and  seasonality  was  used:  this  is  the  Winters  model. 
Some  results  are  shown  in  Table  C.1A. 

Table  C.2  is  developed  from  Tables  C.l  and  C.1A.  It  lists  the  feasible 
exponential  smoothing  models,  based  on  one-step-ahead  forecast  errors.  Table 
C.2  should  serve  as  only  a  guide. 

It  may  be  of  interest  to  note  that  the  smoothing  constants  were  found  by 
means  of  a  computer  search.  In  effect,  the  smoothing  constant  was  chosen  so  as 
to  minimize  the  mean  absolute  deviation  (MAD)  of  the  historical  data. 

It  may  also  be  of  interest  to  remark  that  for  most  of  the  tonnage  series, 
the  exponential  smoothing  models  are  robust  to  "data  updates."  For  example, 
the  percentage  forecasting  error  has  not  changed  considerably  (except  for 
CONSU — Series  //17)  when  the  actuals  for  July  or  3rd  Otr.  1978  were  updated. 


RESULTS  OF  EXPONENTIAL  SMOOTHING — BROWN  1 S  MODEL 
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RESULTS  OF  EXPONENTIAL  SMOOTHING  (WINTERS’  MODEL 


TABLE  C.2 


FEASIBLE  EXPONENTIAL 

SMOOTHING 

MODELS 

Brown 1 s 

Winter's 

Smoothing 

One-Step-Ahead 

One-Step-Ahead 

Series 

Model 

Constant 

Forecast  Error 

Forecast  Error 

PAC/AIR 

Single 

0.60 

0. 72 

(Chopped) 

Linear 

0.30 

1.3 

Quadratic 

0.06 

4.8 

ALA/AIR 

Single 

0.24 

-  1.4 

Linear 

0.12 

0.7 

Quadratic 

0.06 

3.9 

NE/AIR 

(Chopped) 

Single 

0.65 

-  3.0 

EUR/AIR 

Single 

0.‘49 

6.0 

Linear 

0.26 

13.5 

Quadratic 

0.10 

-  7.0 

SOUTH/AIR 

Single 

0.10 

-10.3 

PAC/SEA 

Single 

0.99 

4.9 

ALA/ SEA 

Single 

0.42 

16.5 

Linear 

0.36 

11.2 

Quadratic 

0.25 

6.7 

EUR/SEA 

Single 

0.42 

-13.2 

22.1 

Linear 

0.15 

-22.2 

Quadratic 

0.09 

-20.6  . 

CONUS 

Single 

0.60 

-19.7 

17.8 

(Chopped) 

Linear 

0.30 

-10.4 

Quadratic 

0.02 

-11.1 

Source.  Table  C.l  &  C.2 
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REGRESSION  ON  FLYING  HOURS 

In  Interim  Report  No.  2,  we  used  all  the  data  that  were  available  and 

regressed  tonnage  on  flying  hours.  The  results  in  terms  of  percentage  error 

were  disappointing,  as  seen  in  Table  D.4.  Consequently,  we  followed  the 

suggestion  that  we  regress  only  the  most  recent  eight  to  ten  observations. 

This  we  did  in  Interim  Report  No.  3.  Table  D.3A  is  a  summary  of  this 

regression  and  Table  D.4A  presents  the  percentage  errors.  For  purposes  of 

comparison  these  tables  are  appended  to  Tables  D.3  and  D.4.  It  is  seen  that 

recent  observations  give  better  results  in  terms  of  percentage  error,  but  that 

2  2 

in  many  of  these  cases  the  R  is  zero  or  near-zero.  A  small  R  says  that 

there  is  no  relationship  between  tonnage  and  flying  hours. 

2 

li/hat  does  it  mean  when  the  R  is  very  small  but  when  the  regression 
equation  developed  gives  good  forecasts?  It  means  that  an  average  value  for 
recent  tonnages  would  by  itself  be  a  good  forecast;  that  in  these  circumstances 
the  knowledge  of  flying  hours  does  not  help. 


\ 


TABLE  D.l 


CODE  USED  IN  REGRESSION  ON  FLYING  HOURS 

MILITARY  AIRLIFT  COMMAND  TONNAGE  MONTHLY  DATA 


PAC  Cl 

ALA  C2 

NE  C3 

EUR  C4 

SOUTH  C5 


MILITARY  AIRLIFT  COMMAND  TONNAGE  AGGREGATED  QUARTERLY  DATA 


PAC/AIR  C17 
ALA/AIR  C18 
NE/AIR  C19 
EUR/AIR  C20 
SOUTH /AIR  C21 


MILITARY  SEALIFT  COMMAND  TONNAGE  QUARTERLY  DATA 


PAC/SEA  C6 

ALA/SEA  C7 

NE/SEA  C8 

EUR/SEA  C9 

SOUTH/SEA  CIO 


CONTINENTAL  US  (CONUS)  TONS  MONTHLY  DATA 

CONUS  Cll 


OVERSEAS  FLYING  HOURS  (OSFH)  QUARTERLY  DATA 


PAC  C12 
ALA  C13 
NE  C14 
EUR  C15 
SOUTH  C16 


SEA  PLUS  AIRLIFT  COMMANDS  TONNAGE  AGGREGATED  QUARTERLY  DATA 


PAC-AGG  C30 
ALA-AGG  C31 
NE-AGG  C32 
EUR-AGG  C33 
SOUTH-AGG  C34 


TABLE  D. 2 


QUARTERLY  DATA:  REPRESSING  TONNAGE  ON  OVERSEAS  FLYING  HOURS  (OSFH) — 

ALL  DATA 


Tonnage 

Overseas 
Flying  Hrs. 

Visual 

Regression 

Models 

Series 

Series 

Start  Date 

End  Date 

N 

Identification 

to  Use 

PAC/AIR 
C17  vs. 

PAC/OSFH 

C12 

1st  Qtr.  '73 

3rd  Qtr.  '78 

24 

Linear  Trend 

Straight  Line 

ALA/AIR 
C18  vs. 

ALA/OSFH 

C13 

It 

It 

II 

No  Discernible 
Relationship 

NE/AIR 

C19  vs. 

NE/OSFH 

C14 

II 

II 

23 

Linear  Trend  Straight  Line 

(Delete  first  quarter  outlier) 

EUR/AIR 
C20  vs. 

EUR/OSFH 

C15 

II 

•  1 

24 

Linear  Trend 
(?) 

Straight  Line 

SOt’TH/AIR 
C21  vs. 

SOUTH /OSFH 
C16 

It 

II 

II 

No  Discernible 
Relationship 

— 

PAC/SEA 

C6  vs. 

PAC/OSFH 

C12 

II 

It 

II 

Quadratic 

a.  Straight  Line 

b.  Parabola 

ALA/SEA 

C7  vs. 

ALA/OSFH 

C13 

II 

II 

II 

Linear  Trend 

Straight  Line 

NE/SEA 

C8  vs. 

NE/OSFH 

C14 

II 

II 

II 

Linear  Trend 
(?) 

Straight  Line 

EUR/ SEA 

C9  vs . 

EUR/OSFH 

C15 

II 

II 

11 

No  Discernible 
Relationship 

— 

SOUTH /SEA 
CIO  vs. 

SOUTH/OSFH 

C16 

II 

II 

II 

Linear  Trend 

Straight  Line 

PAC-AGG 
C30  vs. 

PAC/OSFH 

C12 

II 

II 

II 

Quadratic 

a.  Straight  Line 

b.  Parabola 

AI.A-AGG 
C31  vs. 

ALA/OSFH 

C13 

II 

II 

II 

Linear  Trend 

Straight  Line 

NE-ACG 

C32  vs. 

NE/OSFH 

C14 

II 

II 

23 

Linear  Trend  Straight  Line 

(Delete  first  quarter  outlier) 

EUR-AGG 
C33  vs. 

EUR/OSFH 

C15 

II 

II 

24 

No  Discernible 
Relationship 

— 

SOUTH-AGG  SOUTH /OSFH 

C34  vs.  C16 

II 

II 

II 

Linear  Trend 

Straight  Line 

CONUS-ACC  Worldwide  Flying  " 

Hours  (WFH) 

It 

II 

No  Discernible 
Relationship 

27 


TABLE  D. 3 


QUARTERLY  DATA:  REGRESSION  EQUATIONS,  TONNAGE  VS.  OVERSEAS  FLYING  HRS.— 

Tonnage, 

Y 

Overseas 
Flying  Hrs. , 

X 

ALL  DATA 

Adjusted 

R2 

Regression  Equation 

Sy 

PAC/AIR 

PAC/OSFH 

94.3% 

Y  =  4055  +  0.0905X 

1,882 

C17  vs. 

C12 

(Recommended) 

ALA/AIR 

ALA/OSFH 

No  discernible  relationship 

C18  vs. 

C13 

NE/AIR 

NE/OSFH 

50.2 

Y  =  96.3  +  0.704X 

263 

C19  vs. 

C14 

(Recommended) 

EUR/AIR 

EUR/OSFH 

3.2 

Y  =  6415  +  0.0323X 

815 

C20  vs. 

C15 

(Not  recommended) 

SOUTH/AIR 

SOUTH/OSFH 

No  discernible  relationship 

C21  vs. 

C16 

PAC/SEA 

PAC/OSFH 

61.4 

Y  =  56718  +  0.300X 

19,917 

C6  vs. 

C12 

(Recommended;  second  order  polynomial 

not  statistically  significant) 

ALA/SEA 

ALA/OSFH 

40.9 

Y  =  -2420  +  0.740X 

1,343 

C7  vs. 

C13 

(Recommended) 

NE/SEA 

NE/OSFH 

3.6 

Y  =  -905  +  2.28X 

3,719 

C8  vs . 

C14 

(Not  recommended) 

EUR/ SEA 

EUR/OSFH 

No  discernible  relationship 

C9  vs. 

C16 

SOUTH/SEA 

SOUTH/OSFH 

39.0 

Y  =  2427  +  - . 873X 

1,426 

CIO  vs. 

C16 

(Recommended) 

PAC/AGG 

PAC/OSFH 

71.1 

Y  =  60773  +  0.391X 

20,949 

C30  vs. 

C12 

(Recommended;  second  order  polynomial 

not  statistically  significant) 

ALA/AGG 

ALA/OSFH 

39.1 

Y  =  -460  +  0.729X 

1,370 

C31  vs. 

C12 

(Recommended) 

NE-AGG 

NE/OSFH 

(No  statistically  significant  linear 

relationship 

C32  vs. 

C14 

found) 

..  TABLE  D.3  (Continued) 

S-  ARTERLY  DATA: - REGRESSION  EQUATIONS.  TONNAGE  VS.  OVERSEAS  FLYING  HRS 

ALL  DATA  “  “  - - 


Overseas 

Tonnage,  Flying  Hrs., 

Y  X 


EUR-ACG  EUR/OSFH 

C33  vs.  C15 


Adjusted 

R2 


Regression  Equation 


No  discernible  relationship 


SOUTH-AGC  SOUT1I/OSFH 

C34  vs.  C16 


43.8%  Y  =  2579  +  0.987X 


1,466 


CONUS-AGG 


Worldwide  Flying 
Hrs.  (WFH) 


No  discernible  relationship 


29 


TABLE  D.3A 

QUARTERLY  DATA:  REGRESSION  EQUATIONS,  TONNAGE  VS.  OVERSEAS 
FLYING  HOURS— RECENT  OBSERVATIONS 


No.  of 


Tonnage 

Y 

Observations 

N 

_r 

£ 

Regression  Equation 

s 

y.x 

PAC/AIR 

C17  vs.  C12 

8 

-0.184 

0.03 

Y  = 

6,672  -  0.0293X 

323 

ALA/AIR 

C18  vs.  C13 

8 

-0.295 

0.09 

Y  = 

2,310  -  0.0707X 

181 

NE/AIR 

C19  vs.  C14 

8 

0.631 

0.40 

Y  = 

207  +  0.489X 

93 

EUR/AIR 

C20  vs.  C15 

10 

0.649 

0.42 

Y  = 

3,963  +  0.0723X 

695 

SOUTH/AIR 
C21  vs.  C16 

8 

0.117 

0.01 

Y  = 

294  +  0.0123X 

53 

PAC/SEA 

C6  vs.  C12 

7 

-0.645 

0.41 

Y  = 

259,131  -  6.15X 

1,288 

ALA/ SEA 

C7  vs.  C13 

10 

-0.102 

0.01 

Y  = 

4,852  -  0.144X 

1,092 

NF./SEA 

CS  vs.  C14 

8 

0 

0 

Y  = 

-65  +  1.76X 

2,911 

EUR/SEA 

09  vs.  C15 

8 

-0 . 504 

0.25 

Y  = 

98,730  -  0.798X 

10,063 

SOUTH/SEA 
CIO  vs.  C16 

8 

-0.476 

0.23 

Y  = 

6,814  -  1.55X 

1,457 

PAC/AGG 

C30  vs.  C12 

7 

-0.641 

0.41 

Y  = 

265,000  -  6.15X 

10,392 

ALA/AGG 

C31  vs.  C13 

8 

-0.001 

0 

Y  = 

5,217  -  0.0022X 

1,378 

NE/AGG 

C32  vs.  C14 

8 

0.116 

0.01 

Y  = 

142  +  2.25X 

2,970 

EUR/AGG 

C33  vs.  C15 

8 

-0.455 

0.21 

Y  = 

100,000  -  0.688X 

9,912 

SOUTH/AGG 
C34  vs.  C16 

8 

-0.476 

0.23 

Y  = 

7,108  -  1.54X 

1,444 

CONUS/ACG 

8 

0.226 

0.05 

Y  = 

4,031  +  0.002 IX 

700 

vs.  WFH 


Note 

2 

The  R  is  not  corrected  for  degrees  of  freedom 


L  ■ 


r 


TABLE  D. 4 

PREDICT  ION  OF  1978  3RD  QTR  TONNACE  USING  OVERSEAS  FLYING  HOURS  (OSFH)— AI.l.  DATA 


Quarterly 
Sor ios 

K,'il 

resslon  Equation 

3rd  Qtr  ’78 
OSFH 
(X) 

Predicted 
3rd  Qtr  ’78 
Tonnage  (Y’) 

Actual 

3rd  Qtr  ’78 
Tonnage 

t 

Error 

VAC/AIR 

V  = 

4055  +  0.0905X 

34,865 

7120 

0.943 

6033 

-18.0 

NL/AIR 

Y  = 

96.3  +  0.704X 

1,476 

1135 

0.502 

994 

-14.2 

LCR/AIR 

V  = 

6415  +  0.0323X 

64,399 

8495 

0.032 

9664 

2.0 

P AC /SEA 

Y  = 

56,718  +  0.300X 

34,865 

67178  (? ) 

0.614 

35054 

-91.6 

A LA /SEA 

Y  = 

-2420  +  0.740X 

9,365 

4510 

0.409 

3823 

-18.0 

N'K/S'EA 

Y  = 

-905  +  2.28X 

1,476 

2460 

0.036 

318 

-- 

l-c,'  i'H/SEA 

Y  = 

2427  +  0.873X 

1,363 

3617  (?  ) 

0.390 

6927 

47.8 

PAL/ ACC, 

Y  = 

60773  +  0.391X 

34,865 

74405 

0.711 

41087 

— 

■  C/AGC 

Y  = 

2579  +  0.985X 

1,363 

3921 

0.438 

2  934 

-33.6 

f 


31 


TABLE  D.4A 

PREDICTION  OF  1978  3RD  QUARTER  TONNAGE  USING  OVERSEAS 
FLYING  HOURS  (OSFH)— RECENT  OBSERVATIONS 


Quarterly 

Series 

Regression  Equation 

3rd  Qtr .  '78 
OSFH 
(X) 

Predicted 
3rd  Qtr.  '78 
Tonnage  (Y) 

r! 

Actual 

3rd  Qtr.  '78 
Tonnage 

1 

Error 

PAC/AIR 

Y 

= 

6,672  -  0.0293X 

34,865 

5,650 

0.03 

6,033 

6.3 

ALA/AIR 

Y 

= 

2,310  -  0.0707X 

9,365 

1,648 

0.09 

1,807 

8.8 

NE/AIR 

Y 

= 

207  +  0.489X 

1,476 

928 

0.40 

994 

6.6 

EUR/AIR 

Y 

= 

3,963  +  0.0723X 

64,339 

8,614 

0.42 

8,664 

0.6 

SOUTH/AIR 

Y 

= 

294  +  0.0123X 

1,363 

311 

0.01 

278 

-11.9 

PAC/SEA 

Y 

= 

259,131  -  6.15X 

34,865 

44,693 

0.41 

35,054 

-27.5 

ALA/ SEA 

Y 

= 

4,852  -  0.144X 

9,365 

3,505 

0.01 

3,823 

8.3 

NE/SEA 

Y 

= 

-65  +  1.76X 

1,476 

2,533 

0 

318 

-- 

EUR/SEA 

Y 

= 

98,730  -  0.798X 

64,339 

47,387 

0 

40,047 

-18.3 

SOUTH/ SEA 

Y 

= 

6,814  -  1.55X 

1,363 

4,701 

0.23 

7,231 

35.0 

PAC/AGG 

Y 

= 

265,000  -  6.15X 

34,865 

50,580 

0.41 

41,087 

-23.1 

ALA/ AGG 

Y 

= 

5,217  -  0.0022X 

9,365 

5,196 

0 

5,630 

7.7 

NE/AGG 

Y 

= 

142  +  2.25X 

1,476 

3,463 

0.01 

1,303 

-165.8 

EUR/AGG 

Y 

- 

100,000  -  0.688X 

64,399 

55,693 

0.21 

48,863 

-14.0 

SOUTH/ AGG 

Y 

= 

7,108  -  1.54X 

1,363 

5,009 

0.23 

2,934 

-70.7 

CONUS/AGG 
Vs.  WFH 

Y 

— 

4,031  +  0.0021X 

709,030 

(WFH) 

5,519 

0.05 

4,857 

-13.6 

Notes 


WFH  2;  Worldwide  Flying  Hours 
2 

The  R  is  not  corrected  for  degrees  of  freedom 
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ANNEX  E 

BOX-JENKINS  MODELS 

Of  the  Box-Jenkins  models  attempted  the  most  powerful  group  was  the  first- 
order  autogressive  model  (1,0,0).  The  results  for  this  model  are  given  in 
Table  E.l.  The  one-step-ahead  percentage  error  for  this  model  appears  reasonable, 
in  most  cases.  The  Box-Jenkins  (1,0,0)  model  was  not  applicable  to 
series  // 8,  // 9,  and  #10  (NE/SEA,  EUR/SEA,  SOUTH/SEA)  because  of  some  obvious 
seasonal  patterns  in  these.  See  Annex  F  for  an  explication  on  the  seasonal 
patterns  in  some  of  the  Sealift  series. 

We  also  attempted  some  other  Box-Jenkins  models  in  the  hope  that  we  would 
improve  upon  the  (1,0,0)  model.  In  some  cases  we  could.  In  series  it 6  (PAC/SEA), 
for  example,  we  reduced  the  one-step-ahead  error  from -26. 3  percent  for  the 
(1,0,0)  model  to  -16.7  percent  for  the  (4,1,1)  model. 


OTHER  BOX-JENKINS  MODELS 


Constant  term  included 
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ANNEX  F 

TIME  SERIES  DECOMPOSITION  FOR 
NE/SEALIFT  (Series  #8)  and  SOUTH/ SEALIFT  (Series  010) 

The  quarterly  data  of  the  Sealift  series  exhibits  special  problems: 

1)  Some  of  these  series  appear  to  be  seasonal — witness  the  spikes  every  first 
quarter  for  series  #8  (NE/Sealif t) ;  and  2)  The  change  of  the  fiscal  year  in 
FY77  has  confounded  the  recent  data. 

We  wish  to  use  a  time  series  decomposition  model  on  some  of  these  series, 
particularly  series  #8  (NE/Sealif t)  and  series  010  (SOUTH/Sealif t) .  A  de¬ 
composition  model  easily  yields  a  seasonal  index.  But  since  there  are  two 
first  quarters  in  FY77,  we  have  decided  to  combine  these  into  data  for  just 
one  quarter.  Thus,  the  number  of  the  period  would  be  consistent  with  its 
quarterly  position  in  the  fiscal  year. 

We  now  concentrate  on  NE/Sealift  and  SOUTH/ Sealift.  The  two  series, 
NE/Sealif t  and  SOUTH/Sealif t  have  posed  special  problems:  most  of 
the  conventional  time  series  models  failed  to  forecast  them  accurately.  An 
inspection  of  their  plots  exhibits  seasonality;  but  the  periodicity  of  the 
peaks  and  troughs  becomes  erratic  toward  the  end  of  the  series.  Both  series 
appear  to  lack  a  trend.  (This  is  true  for  SOUTH/Sealif t  if  we  delete  the  data 
for  the  fiscal  year  1973).  A  tentative  way  to  handle  these  two  series  is  to 
calculate  for  them  seasonal  indices  by  means  of  moving  averages.  These  cal¬ 
culations  are  given  in  Tables  F.l  -  F.3  for  NE/Sealift  (series  0 8)  and  in 

Tables  F.4  -  F.6  for  SOUTH/Sealif t  (series  0 10). 

Tables  F.7  and  F.8  give  the  forecasts  and  forecast  error  for  the  third 

quarter  of  1978.  The  forecast  errors  are  high,  and  the  seasonal  indices 
may  have  to  be  adjusted  if  these  errors  persist. 
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TABLE  F.l 


TIME  SERIES  DECOMPOSITION:  NEA/SEALIFT  (Series  #8) 


£tr 

Yr. 

Period 

X 

t 

Four-quarter 
Moving  Average 

4x4  Centered 
Moving  Average 

Percentage  of 
Centered  Moving 
Average 

1, 

'73 

1 

8.715 

- 

- 

2, 

'73 

2 

508 

2,590.75 

- 

- 

3, 

'73 

3 

736 

3,486.50 

— 

- 

4, 

’73 

4 

404 

3,538.75 

3,288.69 

12.28 

1. 

’74 

5 

12,298 

3,538.75 

3,554.94 

345.94 

2, 

’74 

6 

717 

3,655.75 

3,391.69 

21.14 

3, 

’74 

7 

736 

2,833.50 

3,209.38 

22.93 

A, 

’74 

8 

872 

2,809.50 

3,020.69 

28.87 

1, 

’75 

9 

9,009 

2,784.00 

2,767.25 

325.56 

2, 

’75 

10 

621 

2,642.00 

2,822.13 

22.00 

3, 

’75 

11 

634 

3,052.50 

2,986.63 

21.23 

A, 

’75 

12 

304 

3,468.00 

3,142.06 

9.68 

1, 

’76 

13 

10,651 

3,405.75 

3,444.13 

309.25 

2, 

’76 

14 

2,283 

3,850.25 

3,236.25 

70.54 

3, 

’76 

15 

385 

2,310.75 

2,165.67 

17.78 

A, 

’76 

16 

2,082 

1,798.25 

2,898.81 

71.82 

1, 

'll 

17* 

4,493 

3,636.00 

2,908.60 

154.47 

0 

*-  * 

'll 

18 

233 

3,889.25 

3,087.75 

7.55 

3, 

'll 

19 

7,736 

3,027.50 

3,394.25 

227.91 

A, 

'll 

20 

3,095 

3,024.25 

“ 

“ 

5, 

'll 

21 

1,046 

- 

6, 

'll 

22 

220 

- 

- 

Change  in  t  e  fiscal  year  (1st  Qtr.  '75);  combined  tonnage 

Source:  Croxton,  Cowdcn,  and  Klein,  Applied  General  Statistics,  Third  edition, 
Englewood  Cliffs,  N.J.:  Front  ice- lla  1 1 ,  pp.  293-302. 


TABLE  F. 2 

PERCENTAGES  OF  CENTERED  FOUR -QUARTER  MOVING  AVERAGES 


FOR  SECOND  DESTINATION  TONNAGE,  NE/SEALIF?  (Series  f 8) 


21.14 
22.00 
70.54 

7.55  I  227.91 


Source:  Table  F.l 


TABLE  F.3 

ARRAY  OF  PERCENTAGES 


ranks  2  &  3 


index 


First 

Quarter 

345, 

.94 

325, 

.56 

309, 

,25 

154, 

.47 

317, 

.4 

332 

.6 

Third 

Quarter 

Fourth 

Quarter 

227.91 

71.82 

22.93 

28.87 

21.23 

12.28 

17.78 

9.68 

22.1 

20.6 

23.2 

21.6 

Source:  Table  F.2 


TABLE  F.4 


TIME 

SERIES 

DECOMPOSITION: 

SOUTH  SEALIFT 

(Series  //8) 

Qtr . , 

Yr. 

Period 

Five-Quarter 
Moving  Average 

Percentage  of 
Centered  Moving 
Average 

1 

'73 

1 

7,726 

— 

— 

2 

'73 

2 

8,109 

— 

— 

3 

'73 

3 

8,610 

7,239.0 

118.94 

4 

•73 

4 

7,481 

6,588.8 

113.45 

1 

'74 

5 

4,269 

5,917.4 

72.14 

2 

'74 

6 

4,475 

5,023.4 

89.08 

3 

'74 

7 

4,752 

4,162.4 

114.16 

4 

'74 

8 

4,140 

4,008.2 

103.29 

1 

'75 

9 

3,181 

3,832.4 

83.00 

2 

'75 

10 

3,493 

3,485.0 

100.23 

3 

'75 

11 

3,596 

3,251.0 

110.61 

4 

'75 

12 

3,015 

3,549.8 

84.93 

1 

'76 

13 

2,970 

3,849.2 

77.16 

2 

'76 

14 

4,675 

4,363.6 

107.14 

3 

'76 

15 

4,990 

5,677.2 

87.90 

4 

'76 

16 

6,168 

5,817.8 

106.02 

1 

'll 

17* 

9,583 

5,421.8 

176.75 

2 

'll 

18 

3,673 

5,093.6 

72.11 

3 

'll 

19 

2,695 

4,258.8 

63.28 

4 

'll 

20 

3 , 34*1 

3,267.2 

102.82 

5 

'll 

21 

1,994 

— 

— 

6 

'll 

22 

4,575 

— 

— 

*Combined  tonnage — change  of  Federal  fiscal  year 


Source:  Croxton,  Cowdren,  and  Klein,  loc.  clt 


TABLE  F. 5 


PERCENTAGE  OF  CENTERED  FIVE-QUARTER  MOVING  AVERAGES 


FOR  SECOND  DESTINATION  TONNAGE,  SOUTH/SEALIFT  (Series  010) 


Fiscal 

Year 

First 

Quarter 

Second 

Quarter 

Third 

Quarter 

73 

— 

— 

118.94 

74 

72.14 

89.08 

114.16 

75 

83.00 

100.23 

110.61 

76 

77.16 

107.14 

i 

87.90 

77 

176.75 

72.11 

63.28 

Source:  Table  F.^ 


TABLE  F.6 

ARRAY  OF  PERCENTAGES 


ranks  2  &  3 


index 


First 

Quarter 

Second 

Quarter 

Third 

Quarter 

Fourth 

Quarter 

176.75 

107.14 

118.94 

113.54 

83.10 

100.23 

114.16 

106.02 

77.16 

89.08 

110.61 

103.29 

72.14 

72.11 

87.90 

102.82 

— 

— 

63.28 

84.93 

80.1 

94.7 

112.4 

104.7 

81.8 

96.6 

114.7 

106.9 

Source:  Table  F.6 
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TABLE  F. 7 


ACTUAL 

ERROR  IN  TONS  AND  PERCENTAGE 

ERROR 

FOR  NE/SEAL1FT 

(Series  #8) 

Quarter 

Seasonal 

Index 

Forecast 
(x  Multipl led 
by  Index) 

Actual 

3rd  Qtr . 

1978 

Percentage 

Error 

Actual 

Error 

(Tons) 

First 

332.6 

9,801 

— 

— 

— 

Second 

33.6 

666 

— 

— 

— 

Third 

23.2 

684 

318 

-115.1% 

-366 

Fourth 

21.6 

637 

— 

— 

— 

Note:  The 

mean  of  the 

series  used  was 

x  =  2,946.9 

TABLE 

F.8 

ACTUAL 

ERROR  IN  TONS  AND  PERCENTAGE 

ERROR 

~  V 

OR  SOUTIi/SEALIFT 

(Series  // 10) 

Quarter 

Seasona 1 
Index 

Forecast 

Cx  Multiplied 
by  index) 

Actual 

3rd  Qtr. 
1978 

Percentage 

Error 

Actual 

Error 

(Tons) 

First 

81.8 

3,825 

— 

— 

-- 

Second 

96.6 

4,516 

— 

— 

— 

Third 

114.7 

5,363 

7,331 

25.8% 

1 , 868 

Fourth 

106.9 

4,998 

— 

— 

Note:  The  mean  of  the  series  used  was  x  »  4,674.5 


ANNEX  G 


POLYNOMIAL  REGRESSION  VERSUS  TIME 

Tables  G.l  and  G.2  give  the  results  obtained  from  regressing  the  tonnage 

series  versus  time.  Table  6.1  gives  the  results  for  a  first-order  polynomial 

(i.e.,  a  straight  line).  Table  G.2  gives  the  results  for  a  second-order 

polynomial  (i.e.,  a  quadratic).  It  is  not  recommended  to  go  beyond  a  quadratic 

in  developing  forecasting  models  using  time  as  an  independent  variable.  In 

2 

terms  of  goodness  of  fit  as  measured  by  R  ,  PAC/SEA  (Series  #6),  ALA/SEA 
(Series  #7),  EUR/SEA  (Series  #9),  and  SOUTH/SEA  (Series  If  10)  perform  remarkably 
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TABLE  G.l 

STRAICHT  LINE  REGRESSION  VERSUS  TIME  (t) 


No. 

Series 

N 

£ 

Regression  Equation 

Predicted 

Tonnage 

Actual 

Tonnage 

% 

Error 

1 

PAC/AIR 

24 

0 

Y  =  1891  -  0.2t 

1,886 

2,011 

6.2 

ALA/AIR 

72 

0.03 

Y 

= 

645  -  0. 768t 

589 

565 

-  4.2 

3 

NE/A1R 

24 

0.04 

228 

4 

EL’R/AIR 

72 

0 

2,506 

5 

SOUTH/AIR 

34 

0 

Y 

= 

109  -  0.105t 

103 

97 

-  6.2 

6 

PAC/SEA 

23 

O.  76 

Y 

= 

133,000  -  4026t 

36,376 

35,054 

-  3.8 

7 

ALA/SEA 

23 

0.55 

Y 

= 

7080  -  196t 

2,376 

3,823 

37.8 

8 

NE/SEA 

23 

0 

Y 

= 

3960  -  84. 5t 

1,932 

318 

— 

9 

EUR/SEA 

23 

0.66 

Y 

= 

84904  -  1566t 

47,320 

40,047 

-18.2 

10 

SOUTH/SEA 

23 

0.39 

Y 

= 

6672  -  166. 5t 

2,676 

7,231 

63.0 

11 

CONUS 

45 

0.03 

Y 

= 

1926  -  1. 24t 

1,869 

1,475 

-26.7 
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TA8LE  G.2 

QUADRATIC  REGRESSION  VERSUS  TIME  (t) 

-  Predicted  Actual  % 


No. 

Series 

N 

Regression  Equation 

Tonnage 

Tonnage 

Error 

1 

PAC/AIR 

24 

0.07 

Y  =•  2010  -  27. 8t 
+  l.lt2 

2,003 

2,011 

0.4 

2 

ALA/AIR 

72 

0.23 

Y  =  542  +  ?59t 
-  0.114t2 

489 

565 

13.5 

3 

NE/AIR 

24 

0.01 

228 

4 

EUR/AXR 

72 

0.13 

2,506 

5 

SOUTH/AIR 

54 

0.02 

Y  =  101  +  0. 732t 
-  0.0152t2 

67 

97 

30.9 

6 

PAC/SEA 

23 

0.77 

Y  =  135000  -  4 512t 
+  20. 3t2 

38,405 

35,054 

-  9.6 

7 

ALA/SEA 

23 

0.63 

Y  =*  8414  -  517t 
+  13.4tZ 

3,724 

3,823 

2.6 

8 

NE/SEA 

23 

0.02 

Y  =  4002  -  94. 5t 
+  0 . 4 1 7 1 2 

1,974 

318 

— 

9 

EUR/SEA 

23 

0.70 

Y  =  78486  -  25.lt 
-  64. 2t2 

40,904 

40,047 

-  2.1 

10 

SOUTH/ SEA 

23 

0.54 

Y  =  8434  -  589t 
+  17.62tZ 

4,436 

7,231 

38.7 

11 

CONUS 

45 

0 

Y  =  1900  -  2 . 18t 
-  0.0743t2 

1,843 

1,475 

-24.9 
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ANNEX  H 

THE  WIENER-KOLMOCOROV  METHOD 

The  Wiener-Kolmogorov  (WK)  method  is  a  spectral  autoregressive  method.* 

It  is  used  here  because  in  theory  it  is  superior  to  other  methods  such  as 
autoregressive  regression.  One  of  its  disadvantages,  however,  is  that  it 
requires  a  long  series  for  parameter  estimation.  (The  same  is  true  of  the  Box- 
Jenkins  model;  but  one  can  always  choose  the  latter  by  guessing  or  by  trial  and 
error) . 


*See  Chan,  K.  H.  and  J.  C.  Havya ,  "Spectral  Analysis  in  Business  Fore¬ 
casting:  The  Wiener-Kolmogorov  Method,"  Decision  Sciences,  Vol.  9,  No.  4 

(October  1978),  pp.  700-711. 


RESULTS  OBTAINED  USING  THE  WIENKR-KOLMOGOROV 

.  -  method 
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Note:  The  numbers  in  parentheses  refer,  respectively,  to  forecasts  and  forecast  errors  using  only  the  mean  and 
trend  components  (i.e.,  the  first  two  terms  of  the  formula).  This  we  shall  call  the  "truncated  VK  model. 


